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Abstract 

The distortion of images of faint background galaxies by (weak) gravitational lensing can 
be used to measure the mass distribution of the deflector. Reconstruction methods for the 
mass proflle of lensing clusters have been developed and successfully tested. Alternatively, 
the image distortions can be used to deflne a weighted mean of the mass inside a circular 
aperture, as was flrst suggested by Kaiser. This "aperture mass" has the advantage that 
strict error bars can be obtained from the data itself, and that a strict lower limit of the 
lens mass inside a circle can be obtained. 

The aperture mass can thus be used to detect dark matter concentrations. Keeping 
in mind that wide-fleld cameras will become increasingly available, this method can be 
used to search for mass concentrations on wide-fleld images. To do this, the aperture 
mass measure is generalized to account for different weighting functions. For each such 
weighting function, a signal-to-noise ratio can be calculated. For an assumed mass proflle 
of the density concentrations, the weighting function can be chosen such as to maximize 
the resulting signal-to-noise ratio. Assuming that dark matter halos can be approximated 
by an isothermal proflle over a large range of radius, a weighting function is constructed 
which is adapted to this density proflle and which yields a smooth signal-to-noise map. 
Numerical simulations which adopt parameters characteristic of 4-m class telescopes are 
then used to show that dark halos with a velocity dispersion in excess of ~ 600km/s 
can be reliably detected as signiflcant peaks in the signal-to-noise map. The effects of 
seeing and an anisotropic PSF are then investigated and shown to be less important 
than might be feared. It is thus suggested that the method of aperture mass measures 
developed here can be used to obtain a mass-selected sample of dark halos, in contrast 
to flux-selected samples. Shear flelds around high-redshift bright QSOs as detected by 
Fort et al. provide a flrst successful application of this strategy. The simplicity of the 
method allows its routine application to wide-fleld images of suflicient depth and image 
quality. 
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1 Introduction 



Objects in the universe are usually discovered by means of their emitted radiation - the 
brighter the objects are, the easier they are detected. Therefore, all samples of objects 
are flux limited. The inclusion of objects thus depends, among other things, on the 
wavelength of observation. As an example, clusters of galaxies can be selected either by 
X-ray surveys or from optical images. There is no reason to expect that the clusters 
selected by these two methods will share most of their properties; in the former case, 
clusters with a dense intracluster medium will be preferentially included in the samples, 
whereas optical survey are sensitive to clusters which contain many bright galaxies. Given 
that the intracluster gas and the galaxies constitute only a small fraction of the cluster 
mass, which is the quantity cosmologists are most interested in, it is by no means clear 
that flux-limited samples reveal a fair representation of (dark) matter concentrations in 
the universe. 

Gravitational lensing offers the possibility to detect 'unseen' matter concentrations 
just from their gravitational effect on light rays. For example, the quantitative analysis 
of the currently existing surveys for gravitationally-lensed QSOs can rule out a signiflcant 
density of compact lenses in the universe with mass ^ 10^1 Mq (e.g., Kochanek 1993, 
1995; Maoz & Rix 1993). Constraints from QSO variability and spectral characteristics 
have yielded strict upper limits on the cosmological density of compact objects with 
IQ-^M© <M< lO^M© (e.g., Schneider 1993; Dalcanton et al. 1994). Since most groups 
and clusters of galaxies are expected not to be sufficiently compact to yield multiply- 
imaged QSO images, these surveys do not yield signiflcant constraints on their cosmic 
number density. 

On the other hand, the mass distribution in clusters can be investigated with grav- 
itational lensing, either from studying arcs and arclets (see the recent review by Fort 
& Mellier 1994), or from the 'weak' distortion of numerous faint background galaxies 
by the tidal gravitational fleld of the cluster (Tyson, Valdes & Wenk 1990; Kaiser & 
Squires 1993). These techniques have already yielded most interesting results; e.g., the 
core radius of the dark matter in clusters was found to be considerably smaller than 
that normally concluded from X-ray studies (e.g., Miralda-Escude & Babul 1995 and 
references therein), and the mass-to-light ratio in the cluster MS1224 was found to be 
considerably larger than 'typical' estimates in other clusters (Fahlman et al. 1994). This 
object supports the possibility that clusters of rather different mass-to-light ratios exist, 
and that optical survey are (strongly) biased towards the detection of only the brightest 
subset of these clusters. 

These lensing studies of clusters have been performed exclusively on clusters which 
were known before from optical or X-ray surveys. However, Fort et al. (1996) have re- 
cently started to develop a strategy for the detection of mass concentrations which are 
not included in existing samples. Supported by the flndings that there is a statisti- 
cally signiflcant overdensity of (foreground) galaxies and clusters around high-redshift 
bright QSOs (Fugmann 1990; Bartelmann & Schneider 1993, 1994; Rodrigues- Williams 
& Hogan 1994; Benitez & Martinez-Gonzalez 1995; Hutchings 1995; Seitz & Schneider 
1995b), on scales in excess of several arcminutes, they assumed that these apparently 
most luminous QSOs are affected by the magniflcation bias. This magniflcation cannot 
be due to single galaxies, as the relevant scale in this case would be of order arcseconds, 
but must come from more extended mass concentrations. Fort et al. (1996) have tried to 
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discover the presence of (foreground) mass concentrations, i.e., by mapping the (weak) 
shear in the field of the QSO, by investigating the systematic eUipticity component of 
background galaxies. In several of the fields which they imaged, a significant shear field 
was detected, which provides a direct verification of the over density of material near the 
line-of-sight to these QSOs. One might suspect that in this way, one could obtain a 
'mass-selected' sample of groups and/or clusters of galaxies. 

The main drawback of this method is that high-redshift, bright QSOs arc rare, and 
these mass concentrations are selected only along these rare lines-of-sight. On the other 
hand, several wide-field imaging cameras are either installed or planned, and the resulting 
images can be used to systematically search for dark matter concentrations. In this 
paper, a statistical method for the detection of (dark) matter concentrations is described 
which is also based on weak lensing of background galaxies. From the image ellipticities 
of background galaxies, one can derive a shear map, and from that, the surface mass 
density of the lens can be reconstructed (Kaiser & Squires 1993; Schneider & Seitz 1995; 
Seitz & Schneider 1995a, 1996b; Bartelmann et al. 1996). One could then analyze the 
resulting mass maps to search for matter concentrations. However, the relation between 
galaxy ellipticities and the resulting mass map is quite complicated, and no simple test 
for a statistical significance of a concentration has been derived yet (although numerical 
simulations can in principle derive such significance levels). Alternatively, Kaiser (1995) 
pointed out that the mean surface mass density in a circular aperture can be directly 
obtained from the shear, 

C{xi,X2) := k{xi) - R{xi,X2) = 2 J 2 — {lt){x) , (1) 

where k{xi) is the mean dimensionless surface mass density within a circle of radius xi 
around a point, R{xi,X2) is the mean dimensionless surface mass density in an annulus 
of radii xi and X2 > xi around the same point, and (71) (x) is the mean tangential 
component of the shear at distance x from that reference point. Given that k is non- 
negative, the above relation can be used to put a rigorous lower bound on the mass 
contained inside the circle of radius xi around the reference point. 

In this paper, the relation (1) is first generalized to allow for arbitrary weight func- 
tions. If one assumes that the mass concentrations one expects to detect have an approx- 
imately isothermal profile over a significant fraction of their extent, the weight function 
can be optimized for the detection of such objects. The important point to note then is 
that these aperture measures have well-defined and simple statistical properties. This al- 
lows to generate a signal-to-noise map from the ellipticities of the faint galaxies, in which 
sufficiently strong lenses will appear as significant (say, ^ 4-cr) peaks. This method is 
then tested on synthetically generated data sets, and it is shown that it lives up to the ex- 
pectation. In particular, the signal-to-noise map is not corrupted by small-scale deflectors 
(such as individual galaxies), nor by larger-scale deflectors. Perhaps most importantly, 
the signal-to-noise maps are much more weakly affected by seeing and an anisotropic 
point-spread function than one might fear, at least as long as the PSF is fairly stable 
over the 'filter scale'. As will be shown in Sect. 4, 'isothermal' mass concentrations with a 
velocity dispersion in excess of ~ 600km/s can be reliably detected by our new method, 
where this lower limit depends on the seeing and the depth of the observations. We 
therefore conclude that this detection method, which is very easy to implement numer- 
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ically and which is fast, can be used routinely on wide- field images to search for dark 
mass concentrations. 



2 Generalized aperture densitometry 



or 



2.1 Notation 

The dimensionless surface mass density /t(x) as a function of angular position d is related 
to the physical surface mass density S{0) by 

«(x) = -^ , (2a) 

^cr 

with the critical surface mass density 

= —^-^ , (26) 

where the -D's denote the distances to the lens and source, and from the lens to the 
source (we use the same notation as in Schneider, Ehlers &; Falco 1992). k is related to 
the deflection potential V'(x) by a Poisson-like equation, 

V'^il) = 2k , (3a) 

V'(x) = -/ dV«(x') ln(|x-x'|) . (36) 

The shear 7 is a two-component quantity which describes the traceless part of the Hessian 
of V^; in particular, we write the shear as a complex number 7 = 71 -|- 172, with 

71 = ^ (V'll - "022) , 72 = -012 , (4) 

where the indices on ifj denote partial derivatives with respect to Xj. Combining (3) and 
(4), one finds that 

7(x) = -/ d^x' P(x - x')^«(x') , (5) 
where the complex kernel T> is given by 

p(x) = -^i^^i±^^ . (6) 

|x| 

If we write the vector x as a complex number X := xi-\- 1x21 then 

where the asterisk denotes complex conjugation. As was shown by Kaiser & Squires 
(1993), the relation (5) can be inverted to yield 



k(x) = —TZe 



TT 



I d2a;'P*(x-x')7(x') 



(8) 



where kq is an arbitrary constant, and TZe{X) denotes the real part of the complex 
number X. Hence, the surface mass density can be determined from the shear field, up 
to an overall additive constant. 



2.2 Aperture measures 

Consider a point xq and a weight function w{x). We then define the quantity 
m(x„) := / d^x .(X) » (|x - x„|) = / d^x .(x + x„) » (|x|) , 



(9) 



which is the integral of k in a circular aperture around xq, weighted by the function w. 
Inserting (8), one finds 



m(xo) 



-7^e 



TT 



j dV7(x')y"d2a;P*(x + xo-x') «^(|x|) + J x kq w {\^\) 



(10) 

The quantity m does not depend on the undetermined constant kq if one requires that 

da; xw{x) = . (11) 



Then, by changing variables in (10) to y = x' — xq, and inserting the definition (7) of V, 
one obtains 



m(xo) 



-1 



-7^e 



TT 



' d^?/7(y + Xo) / dxxw{x) I 
r2 Jo Jq 



d(/7 



(X - YY 



(12) 



where we have again used capital letters to denote complex number, i.e., Y ~ ?/i+i?/2- The 
final integral in (12) can be evaluated by substituting X = xe"'', so that dip = —idX/X, 
and the integral is transformed into a loop integral. The resulting integrand then has 
two poles, one at X = 0, and one dX X — Y . Depending on whether x < \Y\ or not, 
the second pole lies inside or outside the loop. Application of the residue theorem then 
yields 

dip TT 



f2n 

Jo 



Y2 



[2H(|X| - |y|) - \Y\6{\X\ - \Y\ 



(13) 



lo {X-YY 

where H(x) is the Hcaviside step function, and S{x) is Dirac's delta 'function'. Defining 
the tangential shear 7t(y; xq) at position y relative to the point xq by 



Y 



one obtains from combining (12)-(14) 



7t(y; xo) ^ f 7(y + xo) 



y2 



(14) 



m(xo) = J 



2 7t(y;xo) 



Q(|y|) , 



(15) 



where we have defined 
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Q{x):—2 / dx xw{x) — x'^ w{x) . (16) 
Jo 

Equation (15) states that the w-weighted integral of the surface mass density can be 
expressed as an integral over the tangential shear, weighted by the function Q. This 
result was first derived by Kaiser et al. (1994; see also the recent preprint by Squires and 
Kaiser 1996 for a generalization to non-circular apertures), using Gauss' law. To see the 
usefulness of this relation, consider the case that w{x) is chosen such that 

w{x) =0 for X > R ; (17) 

then, owing to the condition (11), Q{x) — for x > R. Thus, in this case, the quantity 
m(xo) can be expressed as an integral of the tangential shear over a finite circle around 
xq. Furthermore, if w{x) is chosen such that w{x) = const, for x < xi < R, then 
Q{x) — for X < xi. Therefore, m(xo) can then be expressed as an integral of the 
tangential shear over the annulus xi < |x — xo| < -R. The ^-statistics (1) derived by 
Kaiser (1995) is a special case (15), which is obtained by setting w — wk, with (setting 
X2 = R) 

WKix) — 9 for < X < xi , 

TTXf 

WKix) = -T-BT~ — 2T xi<x < R, ^^^^ 

n{R^ — xf) 

WKix) = for X > R . 

The obvious advantage of the ^-statistics (1) is that it can be used to directly infer a lower 
bound on the mass inside a circle of radius xi around a point xq, due to the non- negativity 
of the surface mass density, i.e., k{xi,X2) > 0. The reason why we are interested in the 
generalized aperture mass measure m is that the weight function w can be adjusted such 
that the m-statistics is 'optimized' for the detection of mass over densities, as will be 
shown in the next section. 



3 Signal-to-noise statistics 

The observational estimate of the shear field 7 proceeds by investigating the distortion of 
the images of faint background galaxies. Using the same notation as in Schneider &; Seitz 
(1995) and Schneider (1995), except that here we have defined the sign of 7 differently 
as in the quoted papers, we denote by e^^^ the complex ellipticity of a source, such that 
for an elliptical source with axis ratio r < 1, |e*^^^| = (1 — r)/(l + r), and the phase 
of e^^^ is twice the angle between the major axis and the positive a^i-direction (see also 
Sect. 4.4 below). The (complex) image ellipticity e is obtained from e^^^ and the local 
lens parameters as 

where g = — k) is the (complex) reduced shear} In the case of weak lensing, which 
should mainly be considered here, (7^7. Assuming that the intrinsic orientation of the 



The transformation (19) is valid only if |gf| < 1; if \g\ > 1, a different transformation applies (see 
Seitz &Z, Schneider 1996a) . However, in this paper we shall concentrate mainly on non-critical lenses, 
for which (19) is valid. 
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sources is random, then the expectation value of the average eUipticity of a local sample 
of galaxy images becomes (see Schramm & Kayser 1995; Seitz Sz Schneider 1996a) 



{^)e = 9^1 ■ (20) 

Thus, the eUipticity of a galaxy image is an unbiased estimate of the local shear in the 
case of weak lensing. 

Let n be the number density of galaxy images, and let Xj and ej denote the posi- 
tion vectors and ellipticities of the galaxies. Defining in analogy to (14) the tangential 
eUipticity eti(xo) of galaxy image i relative to the point xq by 

= -'^ (^ff^) . w 

where again complex number representation of the vectors x^ and xq was used on the 
right-hand-side, we can write down a discrete version of (15) as 



m(xo) = - V^^^^g(|x,-xo|) , (22) 

" ^ Xo 



n 

i 



where the sum extends over all galaxy images within the aperture, i.e., within the region 
where Q ^ 0. 

In the case of no lensing, the expectation value of m vanishes, (to)^ = 0, because 
the expectation value of the tangential eUipticity is zero. The dispersion amd of m can 
be calculated by squaring (22) and taking the expectation value, which leads to 

2n ^ |xi-xo| 

where the fact was used that the expectation value of the product of tangential ellipticities 
of different images vanishes, (etietj) = for z 7^ j, since the orientation of galaxy images 
are assumed to be mutually unrelated; is the dispersion of the intrinsic eUipticity 

distribution, ^kl^^ = = "^e; the factor 1/2 is due to the fact that the 

dispersion in one component of the eUipticity is 1/ y/2 times the dispersion of the complex 
eUipticity, {eueu) = Ce/2- Thus, we can define the signal-to-noise ratio S at position xq 

as 

\l |Xi-Xor 

Note that S is independent of the normalization of Q, and thus the normalization of w. 

Up to now the weight function w{x), and thus Q{x), has been undetermined. We 
shall now derive the shape of w{x) which maximizes the signal-to-noise ratio for a given 
mass profile. In order to simplify the notation, we set xq = in the following considera- 
tion. The expectation value of m is 

(™)d = ^ET^^(l^*l) • (25) 

Xj 
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Assuming that w{x) = for x > R, so that Q{x) = for x > R (we shall call R the 
'aperture size', or the 'filter scale'), we can perform an ensemble average by averaging 
(25) over the probability distribution for the galaxy positions, 



n / ("^)d = 27r / dxx (7t) (x) 



Qix) 



X' 



(26) 



where N is the expected number of galaxies within R, N = nyri?^, and (7t) (x) is the 
mean tangential shear on a circle of radius x around the point xq = 0. The indices 'c' and 
'd' stand for continuous and discrete. Accordingly, we can perform a similar ensemble 
average of the dispersion (23), to obtain 



■ N 

n 

.fc=i 



^2 



n 







X"^ 



(27) 



From (26) and (27) one can define the ensemble-averaged signal-to-noise ratio Sc as 



20m /o dx (pit) {x)Q{x) / X 
^J^dxQ^{x)/x^ 



(28) 



According to the Cauchy-Schwarz inequality, 

2 



j dx f{x)g{x)^ < (^j dx f\x)^ (^j dx g\x) 



with equality if and only if f / g — const., one sees by identifying / 
(7t) Q/x that S'^ is maximized if 

Q{x) oc (7t) {x) . 



Q/x^/'^ and fg 



(29) 



The corresponding weight function w{x) can be obtained by differentiation of (16), which 
yields a first-order ordinary differential equation for w; the integration constant is fixed 
by the requirement (11). The result is 



w(a;) 



, (7t) iv) , . 
— (7t) [x) 



(30) 



The result (30) for the weight function w which maximizes the signal-to-noise ratio for a 
given functional form of (7t) {y) can be written is a more suggestive form. For that, we 
use the fact that 

R{y) - {k) {y) = (7t) (y) (31) 

(see Bartelmann 1995). Here, K{y) is the mean mass density within a circle of radius y, 
and (k) (y) is the mean mass density on the circle with radius y. We can actually obtain 
this result, which is trivially valid for axially-symmetric mass distributions, by setting 
w{x) = H(j/ — x)/{7ry^) — 5{x — y)/{27ry) in (9). Then, m becomes the left-hand-side of 
Eq. (31), whereas by inserting the corresponding function Q into (15), the right-hand-side 
of (31) is obtained. Using that 
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2 

k{x) — dx' x' (k) {x') , 

^ Jo 

one finds that 

Jy Jy Jo Jy 

After interclianging tlie order of integration, one finds from (30) after a short calculation 
that 

w{x) = (k) (x) - R{R) . (32) 
Correspondingly, one obtains, either from (32) and (16), or from (29) and (31), that 

Q{x)=x^[K{x)-{K){y)] . (33) 

Whereas this 'optimal' weight function maximizes the signal-to-noise ratio, and should 
therefore be best suited for the purpose of this paper, it is not practical to use it, for 
two reasons: first, the mass profile of density inhomogeneities is not known a priori, and 
thus the optimal weight function cannot be determined. Secondly, the weight function 
is neither smooth nor continuous at x ~ R, and therefore the 5" statistics will not be 
a continuous function of position xq. Since we want to search for significant peaks in 
the ^'-function to identify mass inhomogeneities, this non-continuous weight function will 
cause patterns in this map which are difficult to interprete. We shall come back to this 
point. 

Table 1. For the weight function w as given in (34), the parameter a, b, and c are given for several 
combinations of fi and V2. In addition, lOOG is given, where G = f dyy~^Q^{y) 



1^1 


1^2 


a 


6(z/2 - 1)' 


c 


100 G 


0.05 


0.80 


0.7133 


0.1511 


0.1239 


0.7574 


0.05 


0.85 


0.7818 


0.1450 


0.1203 


0.7667 


0.05 


0.90 


0.8528 


0.1398 


0.1170 


0.7758 


0.05 


0.95 


0.9256 


0.1351 


0.1138 


0.7845 


0.10 


0.80 


0.7148 


0.3596 


0.2555 


2.4449 


0.10 


0.85 


0.7825 


0.3435 


0.2483 


2.4844 


0.10 


0.90 


0.8531 


0.3298 


0.2415 


2.5223 


0.10 


0.95 


0.9257 


0.3180 


0.2351 


2.5588 


0.15 


0.80 


0.7171 


0.6251 


0.3844 


4.8258 


0.15 


0.85 


0.7837 


0.5941 


0.3740 


4.9213 


0.15 


0.90 


0.8535 


0.5682 


0.3642 


5.0125 


0.15 


0.95 


0.9258 


0.5461 


0.3549 


5.0997 


0.20 


0.80 


0.7201 


0.9544 


0.5050 


7.8277 


0.20 


0.85 


0.7851 


0.9022 


0.4922 


8.0142 


0.20 


0.90 


0.8541 


0.8593 


0.4800 


8.1912 


0.20 


0.95 


0.9259 


0.8232 


0.4684 


8.3591 



Nevertheless, the optimal weight function should provide a guide for constructing 
weight functions which give rise to a smooth ^-map. As for the weight function (18) 
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used by Kaiser (1995), we want to have that Q{x) = for a; < xi =: viR - for example, 
if the mass overdensity is associated with bright galaxies, there will be no faint galaxies 
visible in the very inner part. This implies that w{x) = const, for x < UiR. Second, 
since one might expect that many of the mass inhomogeneities will have approximately 
an isothermal profile, corresponding to (k) (x) oc 1/x, the weight function should behave 
like 1/x — const, over a large interval. Third, in order to obtain a smooth function jS'(xo), 
w{x) should go to zero at a; = i?. The following form for w{x) is suggested: 

w{x) = 1 for X e [0, viR] , 

w(x) = I — = — c I for a; e fz^ii?, z^2-Rl , f'iA\ 

l-c\^{x-iy^R)^ + iu,R)^ J Li'2j, (34) 

w{x) = ^(-R ~ x)'^{x — aR) for x e [f2-R, R] , 

where the constants a, h and c are determined from requireing continuity and smoothness 
(i.e., continuous first derivative) at x = z^2-R, and by the normalization constraint (11). 
Whereas these constants can be calculated analytically, the resulting formulae are not 
very enlightening. For several combinations of fi and z^2, they are tabulated in Table 1. 
The function w{x) is plotted for four different combinations of vi and V2 in Fig. 1. 




0.2 0.4 0.6 0.8 1 
x/R 

Fig. 1. The weight function w{x) as defined in (34) as a function of x/R, for {vi,U2) = (0.1,0.8) (sohd 
curve), (0.1,0.9) (dotted curve), (0.2,0.8) (short-dashed curve) and (0.2,0.9) (long-dashed curve) 



The corresponding function Q{x) can be calculated from (16); again, the resulting 
expressions are fairly cumbersome and shall not be reproduced here. For the same values 
of vi and 1^2 as in Fig. 1, the corresponding function Q{x) is plotted in Fig. 2. 

Finally, the special case of a singular isothermal sphere is considered. In this case, 
(7t) = xo/{2x), where 

^^-H^y^ (35) 

is the Einstein radius of the mass distribution. Using this tangential shear distribution 
in (28), one finds 

^ _ ^/^rnxl J^dxQ{x)/x'^ _. ^ (3g^) 



dx Q^{x)/x^ 



10 



' ' I ' ' ' I 




0.2 0.4 0.6 0.8 1 
x/R 

Fig. 2. For the same values of vi and V2 as in Fig. 1, the corresponding function Q{x)/x^ is plotted as 
a function of x/R. The identification of the curves is the same as in Fig. 1 

Noting that the ratio of integrals does not depend on R for the function Q{x) correspond- 
ing to w{x) in (34), we have the remarkable fact that the signal-to-noise for a singular 
isothermal sphere is independent of the aperture size R, but depends only on the expected 
number of galaxies within the Einstein radius xq and the intrinsic ellipticity dispersion 
of the sources. For the weight function (18) defined by Kaiser, Q = for a; < uiR — xi, 
and Q = const, for x G [z^i-R, R]. Also in this case does the ratio s of integrals not depend 
on i?, and one can easily calculate that then s = ^/2{T^^u^JJJT+~l^ . The ratio s is 
plotted in Fig. 3 as a function of ui, for 1^2 = 0.9; in fact, the ratio s depends only weakly 
on z/2, as long as (1 — z/2) ^ 1- As can be seen, reducing ui increases the signal-to-noise 
considerably for the weight function (34); for ui ^ 0.1, the Kaiser weight function yields 
larger signal-to- noise ratios. This is due to the fact that for larger values of 1^1, the 
weight function (34) deviates significantly from the optimal weight function (32). Since 
the choice (34) of w is by no means unique, one might try to construct functions which 
yield a higher signal-to-noise for ones favourite mass distribution. In this paper, we shall 
use w as given by (34), and restrict out attention to small values of ui. 




0.1 0.2 0,3 

^1 

Fig. 3. The ratio s as defined in (36a), as a function of vi. The solid curve is for the weight function w{x) 
as defined in (34), for 1/2 = 0.9, whereas the dashed curve has been calculated for the weight function 
(18) as defined by Kaiser (1995) 
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For reference, we rewrite (36a) in useful units as 



4 Detection of mass concentrations from synthetic data 

In this section the finding of mass concentrations from the shear pattern in faint galaxy 
data is discussed. For this purpose, synthetic data are generated, which consist of a set 
of position vectors of galaxy images and their complex ellipticity e^. From these data, 
a signal-to-noise map is obtained, by calculating the function S'(x) as defined in (24) on 
a grid. Mass concentrations should show up as significant peaks in the 5'-map. 

Galaxies are distributed randomly on a square of side Xf ; assuming a number density 
of n, N = nxj galaxies are placed on this square. For each such galaxy, the intrinsic 
ellipticity e^^^ is drawn from the distribution function 

7rcr^(l — e ^I'^e j 

and the value cXe = 0.2 is taken throughout this paper (see Seitz & Schneider 1996, and 
references therein). For a given lens model, the surface mass density «:(xi) and the shear 
7(xj) is calculated at each galaxy position, and the 'observed' ellipticity ei is calculated 
from (19). By selecting the galaxy positions randomly on the observer's sky, we have 
implicity assumed that the magnification bias can be neglected, i.e., that the number 
counts of galaxies follows approximately the law dlog A^(m)/dm cx 10'*'"^, with 7 pa 0.4. 
For a given filter size i?, the ^'-statistics is calculated on a square of size Xi — 2R, to avoid 
boundary effects. 

4.1 Lens-free S'-statistics 

We first consider the case that no lenses are present, i.e, that = ef'K In this case one 
expects that the ^-statistics (24) has a Gaussian probability distribution of unit width, 

p(5) = -i=e-^V2 . (38) 



To verify this expectation, a simulation was carried out with Xf = 60', a number density 
of n = 30/arcmin^ (this value will be used throughout the rest of this paper), and a grid 
of 1500^ points at which S was calculated. The resulting probability distribution p{S) 
is plotted in Fig. 4, together with the Gaussian (38). As can be seen, the probability 
distribution of S follows the expected Gaussian curve within the numerical errors. 

In addition, a modified 5'-statistics was calculated which uses 5" = m/amc instead 
of 5" = m/amd- Whereas the former is slightly more faster to evaluate, since amc is 
independent of xq. Fig. 4 shows that the probability distribution for S' has somewhat 
broader wings than that of S. The same simulations have been carried out with the 
weight function wm, and the results are basically indistinguishable from those shown in 
Fig. 4. 
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Fig. 4. The probability distribution p{S) is plotted (dashed curve), together with the Gaussian distri- 
bution (38) (solid curve). Also, the distribution function for the S' statistics as defined in the text is 
plotted (dotted curve), which shows somewhat broader wings than the S'-statistics. For these curves, 
R = 0'.5 was used, with i/i =0.1 and 1^2 = 0.9 

To illustrate the spatial distribution of the function 5', Fig. 5 displays contour maps 
of on a square of sidelength 10'; i.e., we have chosen Xf = 10' + 2R. No lensing has been 
assumed. Three different filter scales R are chosen in Fig. 5 for the weight function (34); 
in addition, an /S-map for the filter wk given in (18) has been plotted in the lower right 
panel. Obviously, by increasing the S'-maps become smoother. It should be noted 
that the 5'-map for the Wk weight function is much 'noisier' than the one for the w weight 
function for the same filter size R (the two right-hand panels in Fig. 5). This is because 
the weight function wk, and thus the corresponding function Q, are discontinuous, so 
that single galaxy images moving into the filter can cause sudden jumps of S. In the 
maps with small R, several peaks with S > 3 (or minima with S < —3) can be seen, but 
extrema with jS'l > 4 are very rare, in accord with the result shown in Fig. 4. 

4.2 Single-lens statistics 

Next we study the statistical properties of the signal-to-noise S for a single isolated 
defiector. For this, a mass model is chosen which closely resembles an isothermal sphere 
with a finite core and a finite truncation radius, described by the (axially-symmetric) 
mass distribution 



where Xc is the core radius, xt the truncation radius, and xq is the Einstein angle of the 
corresponding singular isothermal model, given in (35). For the rest of the paper, a fixed 
value of Dds/Ds = 0.7 is assumed. The (tangential) shear corresponding to this mass 
distribution is calculated as 7 = 7t = k — k, where 



For this mass model, the expected signal-to-noise ratio Sc can be calculated from (28); 
for different values of cr„ and the ratio Xq/xq, this expected value is plotted as a function 
of the filter scale R in Fig. 6. In addition, simulations were made by distributing galaxies 
around these lenses with an intrinsic ellipticity distribution (37), then distorting them by 



k{x) 




(39) 



2^yx^+xl 2^x'^+xl 
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Fig. 5. Contour plots of the signal-to-noise ratio S defined in Eq. (24), for a surface number density of 
n = 30/arcmin^, on a square of sidelength 10'. No lensing has been assumed, i.e., the galaxy images are 
distributed randomly with an ellipticity distribution according to (37). In the upper left panel, R = 0'.5, 
in the upper right panel, i? = 1', and in the lower left panel, R = 2' . In these three cases, the weight 
function (34) was used, with = 0.1, 1^2 = 0.9. The lower right panel uses the weight function i^k 
given in (18), for R = 1' and vi — xi/R — 0.3. Contour levels are S — 1, 2, 3, . . . (solid contours) and 
S = —1, —2, —3, . . . (dashed contours) 



the lens using (19), and calculating S according to (24) at the center of the lens, again 
for varying filter scale R. For different realizations of the galaxies, the resulting values 
of S are also plotted in Fig. 6. 

Three different cases are shown in Fig. 6: for x^/xq = 1/3, the lens is critical; it is 
nearly critical for Xc/xq = 1, and a rather weak lens for Xc/xq = 3. The ensemble- 
averaged signal-to-noise Sc as calculated from (28) lies well in the 'middle' of the 
numerically-generated points for the non-critical lens, but significantly underestimates S 
for the critical lenses. This is because the relation g ^ j - see (20) - is no longer valid 
near the center of a critical lens. In the calculation of (28) the true value of 7 is used, 
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Fig. 6. The values of S — see Eq. (24) — for different realizations of the galaxy distribution and different 
filter radii R. The numbers in each panel shows the pair of values {cry,Xc/xo) (with cr„ in km/s) used 
in the lens model (39). The truncation radius is chosen to be xt = 50a:;o. For the weight function Q, 
the values i/i =0.1 and 1^2 = 0.9 were chosen. The solid curve in each panel is the ensemble-averaged 
signal-to- noise Sc as given in (28). The systematic differences between the numerical results (crosses) 
and the analytic expectation are due to the non-linearity of the image distortion - see text 



whereas for calculating S according to (24), et is used. Now, e is an unbiased estimate 
for which is larger in magnitude that 7, so that is biased toward lower values for 
strong or even critical lenses. 

For critical lenses with cr„ = 400 km/s, the mean signal-to-noise is of order 3, as also 
follows from (36b). Hence, these lenses can be detected with a marginal significance, 
given the fact that 5" varies considerably around this expectation value for individual 
realizations. However, if the core radius becomes larger than about xq, these lenses are 
too weak to be reliably discovered individually by the method proposed here. Lenses with 
ay ^ 600 km/s can be discovered reliably if the core size is not too large. In addition 
to the core radius, the filter scale R plays an important role: if this filter scale is not 
considerably larger than the core radius, the ^-statistics becomes very insensitive, since 
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it measures radial gradients. Hence, if the filter has the same size as the core, no gradient 
can be measured. For critical lenses, even the smallest filter scale considered in Fig. 6 is 
sufficient, whereas for weak lenses with = Sxq, the filter scale R must be larger than 
2'.5 to detect reliably a lens with — 800km/s. For such large core size, lenses with 
ay = 600km/s are not reliably detected, but those with Xc = xq are ii R ^ 1'.5. Note 
that the width of the distribution of points in all panels of Fig. 6 is of order unity. 

4.3 S'-statistics for an ensemble of lenses 

Next we want to investigate whether the S'-statistics is significantly degraded by not 
considering isolated lenses, but by studying an ensemble of lenses. Two kinds of effects 
may affect the S'-statistics: small-scale mass concentrations which can locally dominate 
the shear field and which should increase the noise in the S'-statistics, and larger-scale 
mass perturbations. The latter should not provide a problem for the detection of mass 
concentrations as long as their contribution can be described as a linear function across 
the size of the filter (i.e., where the weight function Q{x) does not vanish). In order to 
study these effect, we have made simulations with rather extreme perturbations. Mo- 
tivated by Fig. 6, it should be investigated whether lenses with the density profile (39) 
and with ay = 600km/s and Xc = Xq can be significantly detected; those will be called 
'main lenses' in the following. They were distributed randomly with a number density of 
16/deg^. In addition, small-scale lenses were added with ay = 220km/s and Xc = xq, and 
a density which is 100 times higher than that of the main lenses. To simulate larger-scale 
perturbations, lenses with ay = 800km/s and Xc = 5xo were distributed with twice the 
density of the main lenses. Judging from Fig. 6, these 'perturbing' mass distributions 
should not be reliably identifyable in the S-map. 

As can be seen from Fig. 7, several peaks occur with S ^ 5, and all of them corre- 
spond to one of the 'main lenses'; in fact, there is no single peak with S > 5 which is not 
located at a 'main lens' position. Some of the large-scale perturbing mass distributions, 
indicated by triangles in Fig. 7, produce local peaks with 3 ^ S ^ 5, but these peaks 
cannot be confused with those corresponding to the 'main defiectors'. From simulations 
of a larger angular field, the statistics of peaks - as shown in Fig. 8 - clearly indicates that 
all 'main lenses' are identified as S > 5 peaks, and none of the peaks with S > 5 is a 'false 
detection'. We can thus conclude that the aperture method for the detection of matter 
concentrations does work not only for isolated axially-symmetric mass distributions, but 
also if one accounts for large- and small-scale 'perturbations' of the lens mapping. 

4.4 Effects of a point-spread function 

One of the most difficult problems for weak lensing studies from ground-based observa- 
tions is the proper accounting for a point-spread function (PSF). The observed image 
ellipticity e*^"'^^-' deviates from the true image ellipticity e, depending on the size of the 
PSF and the size of the image. Small elliptical images are transformed into nearly round 
observed images by an isotropic PSF. In addition, anisotropics of the PSF can intro- 
duce spurious shear signals in the data. For these reasons, cluster-mass reconstructions 
from ground-based images have to correct for PSF effects (e.g.. Bonnet & Mellier 1995; 
Kaiser, Squires & Broadhurst 1995). Depending on the seeing conditions, rather large 
correction factors have to be used to translate the 'observed shear' to the 'true shear'. 
In this subsection, it will be demonstrated that PSF effects arc much less important for 
the detection of dark matter concentrations if the method described before is applied. 
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Fig. 7. A map of the 5-statistics, with fiher scale R = 2' , for the lens distribution as described in the 
text. Contours are 5 = 1,2, 3, . . ., crosses denote the positions of the 'main' lenses with ay = 600km/s 
and Xc = XQ, and triangles the position of the 'perturbing' lenses with av = 800km/s and Xc = 5xq 



In order to investigate the effects of a PSF on tlie observed image ellipticity, we note 
that the complex ellipticity is defined in terms of the second moment tensor Qij of the 
image brightness distribution as 



<5ll - <522 + 2iQi2 
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where 



Qij — 



J d^x /(x) {xi — Xi){xj — Xj) 



Jd^x /(x) 



(40) 



(41a) 



and X is the 'center' of the image, 
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Fig. 8. Peak statistics: on a field with side = 80', lenses with (t„ = 600km/s and Xc = xq ("main 
lenses") were distributed with a density of 16/deg-^; in addition, to simulate small-scale noise, hundred 
times as many lenses with av = 220km/s and Xc = xq, and for simulating large-scale perturbations, two 
times as many lenses with ctv — 800km/s and Xc — 5a;o were distributed. On the resulting S-map, all 
peaks with S > 3 were identified, and their separation Dmin to the nearest main lens was calculated. 
Plotted are the peaks of the S'-map in the -Dmin-'S'peak plane. Circles corresponds to those peaks which 
are the closest ones to a main lens, and crosses denote the others. Here, a filter scale of ii = 2' was used. 
All peaks with Speak > 5 are situated within R/7 from a main lens; hence, all peaks with Speak ^ 5 
correspond to a mass concentration. Of the 29 main lenses in the field, the lowest one has Speak — 5.36, 
three have Speak G [5,6], nine have Speak G [6,7], 13 have Speak G [7,8], and the other four main lenses 
have Speak > 8. Thus, not only are there no false detections of main lenses for Speak > 5, but also all 
the main lenses are detected with Speak > 5 



_ rd^a;/(x)x 
X = T( \ ■ 



In a similar way, we define the complex ellipticity e*^"^^-* in terms of the second moment 

)(ot 



tensor Q^^^^\ which is defined in terms of the observed brightness profile 



/(°bB)(x) = j d^y /(y) P(x - y) , (42) 

where -P(x) is the PSF, which is assumed to be normalized. From these relations, one 
finds that the relation between Qij and Q^^^^^ reads 



where 



and 



Qf'^ = Qi3+Pi3 , (43) 



Pij = J d^z P(z) {zi - Zi){zj - Zj) , (44a) 
z = j d?z P{z) z . (446) 
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Prom these relations, one can obtain e^°^^^ in dependence of e. A slightly more convenient 
way to find this relation is to consider the ellipticity variable 



Qii - Q22 + 2iQi2 _ 2e 
Q11 + Q22 l + |e| 



X= ^''~^",y^' =T^ , (45) 



and correspondingly for ■)(^(°^^\ Then, in agreement with a similar result derived in 
Bartelmann (1995) and Villumsen (1995), 

X^°'^^ = ^^^^^ , (46) 

where T describes the ratio of the size of the PSF and the size of the image, and x^^^^"* 
accounts for a possible anisotropy of the PSF, 

y ^ -Pii + P22 . (PSF) ^ -Pii - P22 + 2iPi2 
Qii + Q22' ^ Pii + P22 

Hence, (46) expresses the fact that for small T (i.e., size of the PSF small compared to the 
image size), the observed ellipticity is approximately equal to the true ellipticity of the 
image, whereas for large T, the observed ellipticity can be dominated by the ellipticity 
of the PSF. Using the inverse of (45), 

^ (48) 



i + \/i-|xl' 



one can now calculate 6*^°*^®^ as a function of e, for a given PSF and for a given image 
size. We define the size 9 of an image, and the size a of the seeing disk as 

e := VQii + Q22 ; (t'' = Pii + P22 , (49) 

so that T = a^/e^. 

As in the idealized case of negligable seeing, we define in analogy to (22) 



(obs) , 



-, (obs) / \ 

W = -E i'' ^'',2 ^(|xi-xo|) . (50) 



n |Xi - xo 



Whereas m*^"^^-* will deviate significantly from m due to seeing effects, and thus no longer 
is an approximation for the tw-weighted mass density as defined in (9), one can still use 
^(obs) measure local mass over densities; in particular, one can in analogy to (24) 
define the signal-to-noise ratio 

V2 E.#^Q(|x^-xo|) 

s^°''\^o)= iL '"rl.,^ > (51) 
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Fig. 9. The signal-to-noise ratio S^°^^^ as a function of the ratio of the size a of the seeing disk and 
the characteristic size of the galaxy distribution. Here, the size distribution (54a) was assumed. An 
'isothermal' lens model of the form (39) was assumed, with Xc = xq and xt = 50a;o. In the upper panels, 
CTi, = 600km/s, in the lower panels, = 800km/s. The filter size R in the two left panels is i? = 1'.5, 
and ii = 3' in the right panels. The value of S'(°'='^) is plotted as crosses for an isotropic PSF (i.e., 
T^(PSF) _ g-j^ whereas the circles show the values for a PSF with x^^^^^ =0.1. As in the rest of this 
paper, the number density of galaxy images was taken to be 30/sq. arcmin. The solid curve displays the 
ensemble-averaged signal-to-noise ratio, (jn^°^^^'^ /cmc, calculated for an isotropic PSF 



where Ue is the dispersion of the observed eUipticities. Hence, whereas the 'signal' 
uiiohs) ^-jYi ]-,g smaller due to seeing than m, so will be the dispersion, and one might 
expect that the signal-to-noise ratio is not as badly affected as the signal itself. 

For an assumed size distribution of the images, one can define the expectation 
value of m(°^®^ for the center of an axially-symmetric mass distribution (again, we shall 
drop the argument xq in the following) as 
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(obs)\ ^ 1 

/ n 



fc=i 



m 



(obs) 



g(x) 



(52) 



a; 



Note that in this equation, ei^^^ is a function of e^®^ and the reduced shear g{x) through 
(19), i.e., to determine the true image eUipticity e, and of 9 (and the PSF) through the 

transformation (45), (46) and (48) from e to e*^"*^^-*. Because of the dehcate dependencies, 
no further analytic progress can be made, and (52) has to be evaluated numerically. 
Similarly, one can calculate the dispersion of mS'^^^^ in the absence of a lens. 



= I jiopm j d^.w p. (.(')) (4°"')T 



dx Q'^{x) 



(53) 



8 

6 

4 

2 


10 

8 

6 

4 

2 






I I I I I I I 1 1 I I I I I I I I I I I I I I I 



o 



O 



O 



O 



xO O 



I I I I I I I I I I I I I I I I I I I I I I I 





I I I I I I 1 1 I I I I I I I I I I I I I I rh I I I I I I I I I I I I I I I I I 'Pi I I Prl 



01234513345 
u/d. ^/ ^0 



Fig. 10. Same as Fig. 9, except that now the size distribution (54b) instead of (54a) was used 
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For the lens model described by (39), we have calculated the ensemble- averaged 
signal-to-noisc, Sc = (m*^"*^^)) jomci according to (52) and (53) for an isotropic PSF, and 
^(obs) £qj, different realizations of the galaxy distribution, as a function of the ratio ct/^o? 
where 6*0 characterizes the size distribution of the galaxy images. In particular, we have 
chosen the following two distributions: 

..W = ^ , (54.) 

or 

It is obvious that S depends only on the ratio a/9o, since there are no other relevant scales 
in the problem. The results of these calculations are shown in Fig. 9 for the distribution 
function (54a), and in Fig. 10 for the distribution function (54b). 

The symbols in Figs. 9 and 10 show 3^°^^"^ for different realizations of the back- 
ground galaxy population, calculated from (51), whereas the solid curve shows Sc = 
^^(obs)N^ /cr^c: for two values of the velocity dispersion cr^ and two values for the filter 
size R. Crosses denote values of S^°^^^ calculated for an isotropic PSF, whereas the 
circles assume that the PSF is anisotropic, characterized by x^^^^^ = 0.1. As can be 
seen, most clearly from Fig. 10, the distribution of the circles is broader relative to that 
of the crosses; this is due to the fact that for an anisotropic PSF, the contribution of a 
galaxy image to the sum in the numerator of (51) does not only depend on the radial 
separation of the image from the point of interest, but also on its polar coordinate, since 
■;^(PSF) enters et as a term oc cos(2<^). Hence, the polar coordinates of the galaxy images 
provides an additional source of 'Poisson noise' in 

^(obs)^ TheeflFect of a PSF is stronger 
for the size distribution (54b) than for (54a), for fixed which is due to the fact that 
for (54a) the number of 'large' sources is considerably larger than for (54b). As can 
be inferred from the figures, lenses with cr„ = 600km/s can still be reliably detected if 
the seeing is not much larger than the characteristic size of the galaxy distribution. 
However, even for a ^ 56*0 can lenses with cr„ = 800km/s be discovered, independent of 
whether the PSF is isotropic or not. Thus we see that the effect of a PSF is much less 
dramatic than might have been expected. This is mainly due to the fact that not only 
the signal is weakened by a PSF, but also the noise, so that the signal-to-noise ratio is 
much more weakly affected by a PSF than the signal alone. 



5 Discussion 

The generalization (9) and (15) of the aperture densitometry (1), first derived Kaiser et 
al. (1994), no longer has such an intuitive meaning as (1), but this generalization allows 
the definition of a measure which yields the highest signal-to-noise for the detection of 
mass concentrations with an assumed radial profile. In addition, by choosing w appro- 
priately, one can construct a smooth signal-to-noise map from an observed image. The 
weight function (34) used throughout this paper was constructed by assuming the dark 
matter concentrations to have approximately an isothermal profile. Numerical simula- 
tions have then shown that dark matter halos can be detected reliably provided their 
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typical velocity dispersion exceeds about 600km/s (for an assumed density of galaxy im- 
ages of 30/arcmin^, typical of current deep images from 4-m class telescopes); inclusion 
of seeing slightly weakens the lensing signal, and, depending on the ratio of the size of 
the seeing disk and the characteristic image size of faint background galaxies, slightly 
more massive mass concentrations are necessary for a robust detection. Even a signifi- 
cant anisotropy of the PSF does not seriously affect the signal-to-noise ratio as long as 
the PSF is roughly constant across the filter scale R. 

The main advantage of the method proposed here is that the signal-to-noise can be 
directly estimated from the data and thus allows the determination of a significance level 
from the data only. In addition, it can be easily calculated from a list of galaxy images 
which contains positions and image cUipticities. For these reasons, a signal-to-noise map 
can be routinely constructed for all sufficiently deep wide-field images with good imaging 
quality. Since wide-field cameras will become increasingly available, this method will 
allow to perform a survey of mass concentrations selected by mass only (in contrast to 
light)! A mass-selected sample of dark halos is much more closely related to theories 
of structure formation in the universe and can be directly compared with numerical 
simulations, whereas optical or X-ray-selected samples of dark halos can be compared 
with the results of numerical simulations of structure formation only with additional 
assumptions concerning star formation histories and efficiencies, feed-back processes (e.g., 
from supernova explosions), or gas-dynamical processes. Results from lensing in clusters 
like MS 1224 (Fahlman et al. 1994) and C10939-f4713 (Seitz et al. 1996) have yielded quite 
a broad range of mass-to-light ratios, which implies that by 'radiation-selected' samples 
of dark halos those with low mass-to-light ratios will be preferentially included. 

A first attempt to discover dark matter concentrations by weak lensing has been 
successfully made by Fort et al. (1996) by investigating the shear field around high- 
redshift bright QSOs which are suspected to be affected by the magnification bias. The 
method presented in this paper will allow to extend the search for dark matter halos to 
a much larger angular region on the sky. 

Further improvements of the detection method are possible. For example, in Sect. 4.4 
the effects of a PSF were not included in the ^'-statistics explicitly, i.e., not attempt for 
correcting the image ellipticities for a PSF was made. If the PSF is known and stable over 
the filter size i?, one can define in analogy to (50) an aperture measure m^^^'^'^ in which 
the sum over the images is weighted according to the size of the images; large images 
contribute a larger signal than smaller images. The weights can again be optimized to 
yield the largest signal-to-noise ratio. In addition, the detection of weak shear from 
the surface brightness autocorrelation function (Van Waerbeke et al. 1996) provides an 
alternative and completely independent means for measuring the shear. This will be the 
subject of a future publication. 

This work was supported by the "Sonderforschungsbereich 375-95 fur Astro-Teil- 
chenphysik" der Deutschen Forschungsgemeinschaft. 
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